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We consider constraints on a phenomenological dark-matter model consisting of two nearly de- 
generate particle species using observed properties of the Milky Way satellite galaxy population. 
The two parameters of this model, assuming the particle masses are > GeV, are i»k, the recoil 
speed of the daughter particle, and r, the lifetime of the parent particle. The satellite constraint 
that spans the widest range of Vk is the number of satellites that have a mass within 300 pc 
M300 > 5 x 10 6 M©, although constraints based on M300 in the classical dwarfs and the overall 
velocity function are competitive for «k > 50 km s~ . In general, we find that r < 30 Gyr is ruled 
out for 20 km s -1 < Uk < 200 km s _1 , although we find that the limits on r for fixed fk can change 
constraints by a factor of ~ 3 depending on the star-formation histories of the satellites. We advo- 
cate using the distribution of M300 in Milky Way satellites determined by next-generation all-sky 
surveys and follow-up spectroscopy as a probe of dark-matter properties. 

PACS numbers: 



I. INTRODUCTION 

Dark matter is the dominant gravitationally attrac- 
tive component of the Universe [2-g. While there is 
a large set of particle candidates 7H17|. we have no 
idea which of these, if any, constitutes the dark mat- 
ter (although there are some things it clearly cannot be, 
e.g., light neutrinos [l8[ ). However, arguably the most 
popular candidate class is "cold dark matter" (CDM). 
This class of candidate, which includes both axions M, M 
and weakly-interacting massive particles (WIMPs) (19| . 
is called "cold" because it is non-relativistic during ma- 
jor events in the early Universe (freeze-out in the case of 
WIMPs, kinetic decoupling for all). This class is popular 
because it is in many ways the simplest; dark- matter can- 
didates come "for free" in many extensions to the stan- 
dard model of particle physics, are in the early Universe 
at the right abundance in most models, and thereafter 
evolve in a way that is consistent with observations of 
large-scale structure. 

However, CDM is not the only viable dark-matter 
candidate class. A number of observations on smaller 
scales have inspired investigations into dark-matter mod- 
els which reproduce the successes of CDM on large scales 
(corresponding approximately to the scales on which L* 
galaxies are observed and larger) while deviating from 
CDM on the small scales which either lack observations 
or for which observations are difficult to interpret. While 
recent work has focused on finding dark-matter candi- 
dates which may boost the light-lepton density through- 
out the Milky Way to explain unexpectedly high electron 
and positron counts [20|-|22[ , the classic arena in which to 
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play the non-cold-dark-matter game is the distribution 
of dark matter in galaxies. In particular, the mismatch 
between the observed number of satellites of the Milky 
Way and the number of massive subhalos predicted in 
CDM simulations (coined the "missing satellites prob- 
lem" (23j ) has inspired a number of models in which ei- 
ther the phase-space density of dark matter is reduced 
or the small-scale power spectrum is cut off (or both) 
relative to CDM [lalll. 

In this paper, we consider a new set of constraints on a 
class of dark-matter candidate which was originally mo- 
tivated by the missing-satellites problem and the mass 
distribution within dwarf galaxies (25U27| . This class of 
model consists of two nearly degenerate massive dark- 
matter species X and Y, where the masses are related by 
My = Mx(l — e) with e<l. In the simplest scenario, 
X decays to Y and a massless particle which need not be 
a standard-model particle. If e is sufficiently small, the 
Y particle receives a non-relativistic velocity kick — e. 
Unlike most decaying-dark-matter models, we consider 
lifetimes r that are comparable to the age of the Uni- 
verse. Previous work has shown that r > 100 Gyr to 
be consistent with cosmic microwave background obser- 
vations if Vk is relativistic [28j], and r > 30 — 40 Gyr 
for Vk 100 km s _1 in order to remain consistent with 
the observed galaxy-cluster mass function and the galaxy 
mass-concentration relation [2{| [3(| ■ Constraints may be 
tighter if the massless particle belongs to the standard 
model Hi], Hit, but are so far lacking for Vk ^5 100 km s _1 
in the case that the massless particle decays to neutrinos 
or does not interact electromagnetically. 

We reexamine this model in light of its original context, 
the observed population of Milky Way satellite galaxies. 
A flurry of work in the past several years has highlighted 
interesting properties of these galaxies that may shed 
light on both galaxy evolution on the smallest scales and 
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on the nature of dark matter (e.g., [33U48j ]). In particu- 
lar, a number of ultra-faint objects have been found in the 
Sloan Digital Sky Survey (SDSS) H|, which deep pho- 
tometric and spectroscopic follow-up have shown to be 
highly dark-matter-dominated galaxies in the Milky Way 
halo with a nearly constant amount of dark matter within 
the inner ~ kpc. In this work, we show which properties 
of this galaxy population yield robust constraints to the 
decay parameter space (v/., r). We use a hybrid method 
combining semi- analytic dark-matter merger trees and 
star-formation prescriptions with simulations of decay in 
isolated dark-matter halos to determine the properties of 
subhalos and satellite galaxies in a decaying-dark-matter 
cosmology to allow for a comparison with the currently 
known satellite population. For the purposes of this 
work, we define "subhalo" as any distinct dark-matter 
clump within a host dark-matter halo, and "satellite" as 
a subhalo that contains stars. In addition, we will high- 
light how the uncertainty in the evolution of baryons in 
low-mass halo complicates inferences about dark-matter 
properties from Milky Way satellites. We show that the 
resulting constraints on decaying dark matter are com- 
plementary to those obtained in Refs. (29l. l30l]. 

The organization of this paper goes as follows. In Sec. 
HI] we describe the method by which we constrain decay- 
ing dark matter with Milky Way satellites, and in Sec. 
IIII1 we describe the constraints. In that section, we ex- 
plore which properties of the Milky Way satellites may 
robustly constrain the decay parameter space, and show 
how the constraints depend on the star-formation prop- 
erties of dwarf galaxies. In Sec. IIV1 we place our findings 
in the context of other work. 



II. METHODS 

In this section, we describe the key observational prop- 
erties of Milky Way satellites that we will use to constrain 
decay properties, and introduce the method by which we 
calculate the effects of decay on satellite populations of 
Milky Way-like dark-matter halos. We describe the ob- 
servables first because they dictate the requirements for 
the decay simulations. 



A. Observational constraints 

As highlighted in the Introduction, the number of 
known Milky Way satellites has approximately doubled 
(with the exact number a moving target) in the past five 
years due to the advent of sophisticated color-magnitude 
filtering techniques to find low surface brightness galaxies 
in the SDSS [H, H3, Ho[ • These galaxies have a number 
of interesting properties, such as having extremely low 
luminosities (Segue 1 has L s» 300£q ) and being incred- 
ibly dark-matter dominated (again, M/L > 10 within 
the half-light radii) [HI, [52j . The galaxies are pressure- 
supported, and have stellar line-of-sight velocity disper- 



sions ctlos °f only a few km s — , significantly smaller 
on average than the classical Milky Way dwarf galaxies. 
To compare these line-of-sight velocities with theoretical 
predictions for dark-matter halos, ctlos is often converted 
to an estimated u max = max(y / GM(r)/r), the maximum 
circular velocity of the satellit e. Below, we will use the 
relation used by iMadau et al.l [53| , v max — \/3ctlos , to 
compare the u ma x of subsets of the subhalos in our sim- 
ulations to those of real galaxies. Mass-modeling of a 
subset of classical dwarfs indicates that this relation be- 
tween ctlos and v max is reasonable [54[. 

Another interesting, and related, property of the satel- 
lites is the inferred mass enclosed within 300 pc of the 
galaxy centers, M 30 o- Several analyses have indicated 
that this mass is nearly constant among the galaxies 
(M 30 o ~ 10 7 M Q ), even though the luminosities span 
roughly five orders of magnitude [lllilllll. The least lu- 
minous spectroscopically confirmed galaxy, Segue 1, has 
M300 ^ 5 x 10 6 M©. Whether this narrow range of M300 is 
an artifact of the selection function of dwarf galaxies or a 
fundamental limit in star-formation physics is a matter of 
debate [451], but in any case, it means that the minimum 
constraint for a dark-matter model is that it produces 
at least enough subhalos with M300 in the range match- 
ing the observed Milky Way satellites, within a Milky 
Way-mass dark-matter halo. Overshooting the number 
of satellites above the M300 threshold is all right because 
it is generally easier to remove mass from the inner part 
of a halo than it is to add mass. 

The third observed property of these satellites which 
is relevant to this work is that they have stars. Al- 
though an obvious point, it highlights the fact that any 
inference about dark matter from the Milky Way satel- 
lite population depends on galaxy formation and evolu- 
tion. Some authors have estimated the luminosity func- 
tion for the full population of Milky Way satellites (tak- 
ing into account SDSS sky coverage and completeness) 
[13, EH, 56] , and used that to test cold or warm dark mat- 
ter paradigms [571 . |58| . We refrain from using the lumi- 
nosity function to constrain the decay model because it is 
highly sensitive to the poorly understood star-formation 
and feedback processes of small galaxies. Instead, we 
simply consider the fact that the observed satellites nec- 
essarily contain stars. 

In summary, there are three observed or inferred prop- 
erties of the currently known population of Milky Way 
satellite galaxies that we use to constrain the decaying 
dark matter scenario: w max , M300, and the presence of 
stars. In addition, we know how many satellite galaxies 
are currently known, and one can estimate the total num- 
ber of satellite galaxies with properties similar to known 
satellites by taking into account the sky coverage of SDSS 
and its selection function [56| . We discuss the details 
and subtleties of comparing these properties to the sim- 
ulated properties of satellites in a decaying-dark-matter 
cosmology in Sec. Mil 



3 



B. Hybrid decay simulator 

In order to constrain t>k and r, we study the u max , 
M30Q , and stellar properties of subhalos and satellites for 
an ensemble of Milky Way-mass halos. We want to ex- 
plore an ensemble of Milky Way-mass halos for each point 
in Vk—T space in order to get a sense of how likely or un- 
likely it is for the observed satellite population to resem- 
ble the simulated population. In light of this goal and the 
properties of the simulated subhalo and satellite popula- 
tions we use to compare with observations, we construct 
a hybrid technique involving both semi-analytic model- 
ing and TV-body simulations to explore the effects of de- 
cay on subhalos and satellite galaxies of Milky Way-mass 
dark-matter halos. 

The fact that we are interested in long (r 3> 1 Gyr) 
decay times and non-relativistic allows us to use CDM 
initial conditions for our hybrid decay simulator. We use 
merger trees from the Galacticus semi-analytic model, 
which are simulated in ACDM cosmologies, and use pre- 
scriptions within Galacticus to determine the density 
profiles of the dark-matter halos and subhalos in the ab- 
sence of decay [59| . The relevant properties of Galacti- 
cus for this work are summarized in Sec. MB 11 We 
use simulations of decay in isolated dark-matter halos to 
take into account the effects of decay on individual sub- 
halos in the merger trees, as described further in Sec. 
Ill B 21 We populate subhalos with stars according to the 
prescription in Sec. Ill B 31 

This method is much faster to implement and in many 
ways more robust than using cosmological ./V-body sim- 
ulations alone, which at first glance would have been the 
obvious path to take. Cosmological TV-body simulations, 
even zoomed on a particular host halo, have major dis- 
advantages. First, even the highest-resolution simula- 
tions of Milky Way-mass halos can only probe down to 
a few hundred parsecs of the center of the main halo po- 
tential (60l - [63l |. However, while resolution tests exist for 
main halo centers, systematic resolution tests on subha- 
los are lacking. This is a problem if we want to probe the 
mass within subhalos on scales comparable to the demon- 
strated resolution limit on the host halo. However, in 
setting the subhalo properties in the merger tree, we do 
use relations among halo properties that are calibrated 
on large TV-body simulations and in ranges of mass or 
redshift that have not been tested by simulations. Sec- 
ond, each realization of a Milky Way-mass halo at the 
present best resolution takes months of supercomputing 
time to run. Since we want to explore a range of and 
t, and to simulate an ensemble of Milky Way-mass ha- 
los for each set of Dk and t, cosmological simulations are 
highly impractical. 



code 1 [59j | . Only dark sector physics (dark matter merger 
tree construction and subhalo orbital decay via dynami- 
cal friction) was included in these calculations — all bary- 
onic physics in Galacticus was switched off. 

Dark matter merger trees were built using the algo- 
rithm described by Ref. [sij. Standard values of the 
accuracy parameters for this algorithm were used as fol- 
lows: 

[mergerTreeBuildCole2000MergeProbability] = 0.1: 

The maximum probability for a binary merger 
allowed in a single time step. This ensures that the 
probability is kept small, such the the probability 
for multiple mergers within a single time step is 
small; 

[mergerTreeBuildCole2000AccretionLimit] = 0.1: 

The maximum fractional change in mass due to 
sub-resolution accretion allowed in any given time 
step when building the tree. 

Merger trees were resolved down to halos of mass 10 7 M Q . 
Mass accretion below this scale was treated as smooth ac- 
cretion and branches were truncated once they fell below 
this mass. 

Branching probabilities in the merger tree were 
computed using the algorithm of Parkinson et al. [o5| . 
The parameters Go, 71 and 72 of their algorithm were 
set to 0.57, 0.38 and —0.01 respectively as recommended 
by Parkinson et al. [sE|. Additionally, the parame- 
ter [modif iedPressSchechterFirstOrderAccuracy] 
in Galacticus was set to 0.1 to limit the step 
taken in the critical linear theory overdensity 
for collapse in the merger tree building algo- 
rithm. This step was not allowed to exceed 
[modif iedPressSchechterFirstOrder Accuracy] 
times v /2[cr 2 (M 2 /2) - <t 2 (M 2 )], where M 2 is the mass 
of the halo being considered for branching and a(M) is 
the CDM mass variance computed by filtering the power 
spectrum using top-hat spheres. This ensures that the 
first order expansion of the merging rate that is assumed 
in the tree building algorithm is accurate. 

Progenitor halo mass functions from merger trees built 
using this algorithm have been compared with equivalent 
progenitor mass functions measured from the Millennium 
Simulation [6611 an d show excellent agreement with the 
N-body result [III, III . 

The Galacticus code evolves the merging distribu- 
tion of halos forward in time. When one halo merges 
with another, larger halo it becomes a subhalo within 
that larger host. We track only a single level hierarchy of 



1. Merger tree 

Distributions of dark matter subhalo properties at z = 
were computed using the Galacticus semi-analytic 



1 Specifically, vO.9.0, revision 12 of GALACTICUS was 
used. The GALACTICUS model can be downloaded from 
http://sites.google.com/site/galacticusmodel The input 

parameter file used for these calculations is available at 
http: / / www.ctcp.caltech.edu/galacticus / parameters / darkMattcrDecays.xml 
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substructure, i.e. we track only substructures, not sub- 
substructures or deeper levels or the merging hierarchy. 
Therefore, if a merging halo contains its own subhalos 
they will become independent subhalos within the new 
host and will be assigned new merging times (see below) . 

Once a halo becomes a subhalo it is assigned a 
timescale for merging to the center of its host halo due 
to the actions of dynamical friction. We use the dynam- 
ical friction calibration of Jiang et al. [67[ to compute 
dynamical friction timescales, with orbital parameters of 
subhalos selected at random from the cosmological dis- 
tribution found by Benson [68[. Once this timescale for 
merging has elapsed the subhalo is merged into its host 
and no longer exists as an independent entity. 

The properties of the subhalos, including the den- 
sity profile and the mass, are set at accretion time. 
The radius at which the slope of the density profile 
r s = dlog p/dlogr = — 2 relative to the virial radius i? v ir, 
defining the concentration parameter 

c = R vh -/r s (1) 

is set according to Gao et al. [69|. which is calibrated 
using the Millennium Simulation [66j . The subhalo mass 
and virial radius at accretion are set according to the 
virial overdensity criterion of Percival (70j for homoge- 
neous dark-energy CDM cosmologies. 

There are a few caveats to applying these particular 
options for the merger tree. First, the merger tree and, 
for example, the mass-concentration relation as a func- 
tion of redshift have been tested on a small set of simu- 
lations representing a limited set of cosmological param- 
eters and range of halo mass. For example, the redshift- 
dependent concentration and the merger histories have 
been calibrated using the Millennium Simulation, which 
has a relatively high as = 0.9, and for which the dark- 
matter particle mass is ~ 1O 9 M (larger than many of 
the subhalos in our merger trees) [66], |69( . Others have 
found that the mass-concentration relation depends on a 
number of cosmological parameters, in particular as, and 
the redshift-evolution of this relation is still under debate 
[7ll - l73j . Moreover, this relation has not been tested in 
the mass and redshift ranges of some of the subhalos be- 
fore they merge onto a larger halo. The main way we 
mitigate some of these uncertainties is to impose a min- 
imum cut-off in the concentration. Studies have shown 
that c ~ 4 at virialization, so we assign high-redshift ha- 
los c = 4 if the Gao et al. [69| formula indicates c below 
that value [74| . In addition, there is a great deal of scatter 
in the mass-concentration even in the mass and redshift 
range in which it has been studied. The dynamical fric- 
tion formula has been determined using a set of N-body 
cosmological simulations, and also not necessarily on the 
small and early scales that are relevant to this work [67j . 

Finally, the merger tree does not include the effects of 
tidal stripping. Tidal stripping not only destroys some 
subhalos, but also may reduce w max and M300 of the re- 
maining subhalos. The important implication of this fact 
for this work is that our constraints are quite conser- 



vative, as we use the subhalo properties at the time of 
accretion. 

2. Decay simulations 

In order to estimate the effects of decay on the subhalo 
population and the host halos, we use a set of simulations 
of isolated, initially equilibrium CDM halos. The first set 
of simulations we use, with c = 5 and 10, were initially 
presented in Peter et al. [3(|. The halos in these simu- 
lations had initial virial mass M v - a = 10 12 M Q , with the 
virial overdensity defined in Ref. (75j . There were 25 
sets of simulations, scanning decay parameters Wk/w V ir = 
0.077, 0.38, 0.77, 1.54 and 3.85, and r = 0.1, 1, 10, 50, and 
100 Gyr, where u v ir is the virial velocity of the halo. For 
this work, we simulated M v ; r = 10 12 M Q halos with one 
million particles each of mass 10 6 M Q with c = 20 and 
c = 30 for the same sets of Vk/iVir and r in order to span 
a broad range of subhalo concentrations. We simulated 
these halos using a modified version of Gadget-2 (66|. As 
in Peter et al. [30J, we assumed a Navarro-Frenk- White 
(NFW) density profile for the initial matter distribution 
within the virial radius i? v j r of the halo [76|, [77} , 

p<r) = / s ,2 , (2) 
-(1 + -) 

r B \ r s J 

where p s is the scale density, and r s is the scale radius. 
At the level of numerical resolution in our simulations, 
this density profile is nearly indistinguishable from the 
now-preferred Einasto profile (78[. All other aspects of 
the simulation, from the set-up of the initial conditions to 
the simulation and cosmological parameters, is identical 
to that in Peter et al. (30j . 

To apply these simulations to the subhalos of the 
merger trees, as well as to the host halos, we use the 
following strategy. To characterize the subhalos, we use 
r s and p s at the time of accretion, and find the mass and 
concentration of these subhalos at z — assuming that 
the only growth to the subhalos comes from the decrease 
in the global mass density of the Universe as a function 
of time. This sets the z = concentration parameter 
and virial mass, and we estimate the effects of decay us- 
ing these properties. We characterize the isolated decay 
simulations by five parameters: UkAViri which relates the 
recoil speed to the typical speed of particles in the halo 
(and hence, the escape velocity); c, which characterizes 
the depth of the potential well in addition to the typical 
dynamical time scale of particles within the halo (since 
the dynamical time at the half-mass radius depends only 
on c); the decay time scale r; r/R v - lr , the radius at which 
we sample the mass profile of the halo as a function of 
the virial radius; and time t. We divide the mass pro- 
file by the initial virial mass in the simulation, and we 
cut off the mass profile at r = 0.04i? v j r since we have 
found that numerical relaxation becomes a problem for 
smaller radii. For a given subhalo with specific z = 
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properties, we find the subhalo mass profile by interpo- 
lating the simulations in the log-space of Uk/i> v ir, c, r, 
and t. We use the mass profile to determine M300 and 
v maK at z = 0, and iWx at the accretion time (since this 
determines whether or not the subhalo has stars, as will 
be described below). 

We also apply the decay interpolation to the host halos, 
since decays can obviously change the mass of the host 
as well as the subhalos. This is important for the host 
mass cut we employ in Sec. IIII1 

This provides a conservative estimate of the effects of 
decay on the subhalos since the interplay between decay 
and tidal forces are likely to reduce v max and M300 even 
more than if the subhalos were isolated halos. 



3. Populating subhalos with stars 

Since we only care if a subhalo has stars at all, a 
crude model for populating halos with stars is accept- 
able. If a halo is accreted onto a larger halo prior to 
reionization, the halo is allowed to have stars if its max- 
imum circular velocity at accretion u max > 2 km s — , 
which is approximately the lowest threshold at which 
gas may accrete onto halos and cool via collisional in- 
teractions with molecular hydrogen (79T - [8lj |. During and 
after reionization, star formation is suppressed in low- 
mass halos due to various effects related to the strong 
ionizing background radiation 82-86]. To model these 
effects, we allow halos with v max > 38 km s _1 at accre- 
tion onto a larger halo after reionization to host stars; 
those with smaller u max are not allowed to host stars. 
Though this step-function treatment of stellar content in 
halos is crude, it captures the essence of the fact that 
star formation in halos depends on reionization and that 
the stellar content of a subhalo depends on its accretion 
history. This model for populating halos with stars is 
similar to that adopted by Madau et al. [531 ] . 



III. RESULTS 

For each set of («k, r), we selected ~ 100 host ha- 
los with virial masses in the range (0.5 — 2) x I0 12 M Q , 
since this appears to be the plausible range of mass for 
the Milky Way halo (although most estimates favor the 
higher end of this range) [87h97| . We have checked the 
results in this section for dependence on the merger his- 
tory. Specifically, simulations suggest that the Milky 
Way disk could not have withstood a 10:1 merger since 
z ~ 1 [98, 99j]. We found that the results from the entire 
host-halo population were indistinguishable from those 
from the sample of hosts selected to have not had a ma- 
jor merger since z = 1. 

In this study, we consider several subsets of the sub- 
halo populations. The first sample we call the "all nodes" 
sample, as it represents all nodes of the merger tree down 



to the mass resolution of 10 7 M Q . This is the most con- 
servative subhalo sample to use to constrain — r pa- 
rameter space because we ignore dynamical friction and 
we characterize the subhalos by their properties at the 
moment they are accreted onto larger halos. The second 
subset of subhalos we use is the "dynamical friction" sam- 
ple. This contains all subhalos in the merger tree that do 
not sink to the center of larger halos by dynamical fric- 
tion. The survival probability is calculated according to 
Ref. [67[ and is the default setting of Galacticus [59| . 
The mass profile of these subhalos varies continuously 
prior to when they first become a subhalo (in either the 
main branch of the merger tree or a sub-branch). After 
that time, their mass profile remains fixed. Using this 
sample to set constraints is less conservative than using 
the "all nodes" sample because it is not completely clear 
how decays will affect dynamical friction; on one hand, 
decays decrease the mass of both the subhalo and host 
(although the effect on the former is stronger than on the 
latter), reducing the effectiveness of dynamical friction; 
on the other hand, the binding energy of the subhalo de- 
creases, making it easier to shred (although this effect is 
not modeled in Galacticus for CDM halos). 

The differences in the CDM properties of these two 
populations is shown in Fig. [T] We show these proper- 
ties for CDM because they will influence the properties 
of the "all nodes" and "dynamical friction" populations 
once we turn on the decays. What we show in this figure 
is the velocity function, the number of subhalos above 
a maximum circular velocity v max , N(> v max ), for both 
the "all nodes" sample (left) and "dynamical friction" 
sample, with the error bars showing the 25% and 75% 
percentiles of the number of subhalos in the host ha- 
los. The percentiles should be taken as an indication of 
the range of i> max and not interpreted strictly as errors 
because we have chosen an initial host population with 
masses drawn uniformly from the range 1-3 x 10 12 M Q 
instead of choosing the host masses according to a prob- 
ability distribution from existing data on the Milky Way 
host mass. These merger trees were generated for a flat 
ACDM cosmology with cr 8 = 0.9 and fi ro = 0.2725 and 
n s = 0.961. Note that since the velocity function is a 
cumulative function, the percentile bars are highly cor- 
related. Also on the plots are the velocity function of 
known Milky Way satellites, corrected for SDSS sky cov- 
erage (dotted line connecting data points) (62[, and the 
velocity function from two different high-resolution CDM 
simulations of Milky Way-mass dark-matter halos (thin 
solid lines). The upper thin line represents the velocity 
function found in the high-resolution Aquarius A simula- 
tion of a M vir ~ 2 x 10 12 M Q halo with fl m = 0.25, n s = 1, 
erg = 0.9 cosmology [61]. The lower thin line shows the 
velocity function of the Via Lactea II halo, which has a 
similar virial mass as the Aquarius A halo but is sim- 
ulated in a cosmology with O m = 0.238, n s = 0.951, 
erg = 0.74 [53j]. This velocity function has a lower normal- 
ization because the cosmological parameters used tend 
to produce fewer and less dense subhalos than the ones 
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employed by the Aquarius collaboration. The velocity 
functions in the CDM simulations lie below those of our 
merger trees because the simulations necessarily take into 
account tidal stripping. Tidal stripping tends to lower 
Vmax after accretion, and Madau et al. [53| shows that 
the velocity function of v ma x at accretion has approxi- 
mately a factor of five higher normalization than the ve- 
locity function for which stripping is taken into account, 
for u max > 6 km s~ . 

It is apparent that velocity function of the "dynamical 
friction" sample has both an overall normalization lower 
than that of the "all nodes sample" , and that the slope is 
substantially steeper. This is to be expected, as dynam- 
ical friction is more efficient for higher-mass (and hence, 
in general, higher v max ) subhalos. Thus, in the case of 
CDM, we find that the "all nodes" sample contains both 
more and on average more massive subhalos. Moreover, 
given the factor of ~ 5 difference in normalization for 
the velocity function for v max at accretion versus v max at 
z = with tidal stripping included, the "dynamical fric- 
tion" sample is a somewhat better match to the veloc- 
ity function of the two high-resolution CDM simulations, 
taking into account differences in underlying cosmologies 
and host halo masses [6l|, HH . 

Within each of the "all nodes" and "dynamical fric- 
tion" subhalo samples, we determine two different possi- 
ble satellite populations: a z vc — 7 satellite sample con- 
sisting of all subhalos satisfying the star-formation cri- 
terion for 2 rc = 7 given in Sec. Ill B 3[ and a z rc = 11 
satellite sample consisting of subhalos satisfying that cri- 
terion at a reionization redshift z rc = 11. We choose 
the z Ie = 11 reionization redshift because WMAP seven- 
year data suggest tha t the Universe was reionized at 
z Ie = 10.5 ± 1.2 [100j . However, reionization is ulti- 
mately a time- and location-dependent process. Other 
work has indicated that the Milky Way-associat ed re gion 
could have been reionized as late as z rc = 7 0, llOlj . 



A. Number of subhalos and satellites above 

M 30 o = 5 x 1O 6 M 

Our first constraint on and r comes from considering 
the number of simulated satellites above an M300 thresh- 
old M 300 = 5 x 10 6 M Q . This threshold for observed satel- 
lites is apparent in Strigari et al. [4l[. To estimate the 
true number of Milky Way satellites above this thresh- 
old from the observed s atellite population, we extend the 
work of iTollerud et all [HI]. In that work, the authors 
use the sky coverage and selection function of SDSS as 
well as th e subhalo distribution in the Via Lactea simula- 
tion [l02l ] to estimate the number of satellite galaxies out 
to the virial radius (~ 389 kpc for the Via Lactea halo, 
Af vir = 1.8 x 10 12 M Q ). They found that there should be, 
on average, 382 satellites n the Milky Way within that 
virial radius, and a 98% probability that there would be 
at least 292. Since we consider the possibility that the 
Milky Way halo could be up to a factor of four less mas- 



sive (remaining consistent with published estimates of 
the Milky Way halo mass), we must adjust the Tollerud 
et al. [42I ] results to find the minimum number of satellites 
within the Milky Way virial radius. 

We do this by considering the radial distribution of 
the Via Lactea subhalos, rescaling the distribution by 
the virial radius. We find that for possible Milky Way 
halo masses M v - U > 5 x 10 11 M Q , we expect at least ~ 
200 satellites within the Milky Way virial radius. This 
lower bound is somewhat roug h, as we do not do a full 
recalculation of ITollerud et al.l but instead estimate the 
average difference in the parameter /(> r) (defined in 
Tollerud et al. 42]) between a virial mass 5 x 10 n M Q 
and 1.8 x 1O 12 M . As long as the radial distribution of 
subhalos is relatively insensitive to host halo mass, and 
as long as the radial distribution is not a strong function 
of decay, this approach should yield an approximately 
correct estimated minimum number of satellites. 

In order to set conservative limits on — t space, 
we create merger trees with erg = 0.9, which is slightly 
above th e 2— a upper limit from the WMAP seven-year 
data set [lOOt ] . We choose a high cr 8 because structures 
form earlier and have higher densities for high <7g. The 
consequences for our study is that there are higher num- 
bers of subhalos per host halo, and those subhalos have 
higher concentration (and are thus less prone to disrup- 
tion and will tend to have high v max and M300) for high 
<7g. In addition, for low-redshift measurements (includ- 
ing the galaxy power spectrum), decays can masquerade 
as low c 8 . Although high- and low-redshift estimates 
of erg are largely consistent with each other, w e choose a 
high erg for our study to be conservative [ll.[a-[5l [T03l[T04| . 
We find that the number of subhalos and satellites for <rg 
at its mean WMAP seven-year values is up to a factor of 
two less than for the as = 0.9 samples. 

In order to compare to decaying-dark-matter cosmolo- 
gies, in Fig. [21 we show the probability distributions for 
the numbers of subhalos and satellites of the "all nodes" 
and "dynamical friction" samples for CDM, according 
the the merger trees. Again, the distributions are not 
actual probability distributions for the Milky Way, as we 
have not weighted the host halo mass distribution ac- 
cording to the probability distribution of the Milky Way 
mass from observations. The distribution is meant to 
give a sense of the range of possible numbers of subhalos 
and satellites within the Milky Way halo. The top panel 
shows the distribution of subhalos for each sample, and 
the distribution in the number of subhalos satisfying the 
M 30 q > 5x 1O 6 M criterion. We find that there should be 
thousands of subhalos satisfying this criterion for either 
the "all nodes" or "dynamical friction" samples, with a 
factor of ~ 3 more expected in the "all nodes" than "dy- 
namical friction" samples. The middle panel shows the 
number of subhalos satisfying the star-formation crite- 
rion of Sec. Ill B 31 regardless of M3oo- There are propor- 
tionally far fewer satellites for the "dynamical friction" 
samples than the "all nodes" samples because many of 
the high-u max subhalos in the "all nodes" sample have 
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all nodes dynamical friction 




FIG. 1: Velocity function of subhalos. The thick solid lines with error bars represent the CDM velocity functions from the 
merger trees for the "all nodes" (left) and "dynamical friction" (right) samples. The dotted line with error bars show the 
velocity function for subhalos with M300 > 5 x 1O 6 M0. The dotted line connecting data points represents an estimated SDSS 
sky coverage-corrected velocity function for known Milky Way satellites and the thin lines (upper: Aquarius A 6lj; lower: 
Via Lactea II [53}]) represent velocity functions found in high-resolution CDM simulations. 



merged with host via dynamical friction. The lowest 
panel shows the distribution in the number of subhalos 
satisfying both the M300 and star-formation criteria. We 
find little in the way of differences between the middle 
and lowest panels of Fig. [5] because of a mild correlation 
between M300 and v max , and a strong correlation between 
halo formation time and M300, such that high-v max ha- 
los that can form stars after reionization also have high 
M300, and the medium- w max halos that form stars prior 
to reionization form early and thus have high M3oo- 

We illustrate the effects of decay on the subhalo and 
satellite samples for = 30 and 200 km s _1 and r = 20 
and 60 Gyr in Fig. (3J which gives a flavor of what de- 
cay does to the subhalo and satellite populations. When 
Vk = 30 km s~ , the smaller subhalos tend to be dispro- 
portionally affected, since v v - 1Y or u max of the larger sub- 
halos are a bit bigger than v^. Much of the dark matter in 
the small halos is quickly ejected, and the daughter par- 
ticles that remain in the halo are responsible for a fairly 
large (but r-dependent) injection of kinetic energy, which 
tends to drastically reduce the central density. The effect 
is more pronounced for smaller r because the central den- 
sity is extremely sensitive to the decay fraction if many 
or most of the daughter dark-matter particles are ejected 
from the halo after the decay. The reason that the num- 
ber of z re = 7 and z re = 11 satellites is nearly identical 
for = 30 km s -1 is that only the highest-u max halos 
are largely unaffected by the decays. 

For the same r, there are more high-M3oo subhalos 



for higher because the decays start affecting the host 
halos, too. If host halos suffer mass loss due to decay, 
then in order for the host halo to be in our specified 
range, it must have had a higher CDM halo mass. Since 
the number of substructures above a mass threshold is 
correlated with host mass, the hosts at z = that are 
significantly affected by decay and whose 2 = mass lies 
within our target range have more subhalos than if the 
effects of decay were minimal. This effect is noticeable 
if one compares the two upper plots with the two lower 
plots in Fig. H 

One potential issue with our hybrid decay simulator is 
that 300 pc is often lower than the smallest r/R vlI bin in 
the halo mass profile (Sec. Ill B 2p . This means that we 
must extrapolate beyond our simulated data to calculate 
M3oo- In general, this means that we tend to overesti- 
mate M300 , since CDM simulations (as well as our decay 
simulations, in the inner region unaffected by numeri- 
cal relaxation) find that the inner slope of the density 
profile tend s to become less steep the deeper one gets 
in the halo [l05j |. We illustrate the effects of the inner 
radial cut-off of the mass profile on the subhalo and satel- 
lite populations in Fig. [4] with cut-offs of r = 0.01i? V i r , 
0.04i? vir (default), and 0.08i? V ir- For these plots, we set 
«k = 30 km s _1 and r = 20 Gyr. The innermost cut-off, 
r = 0.01i? v i r , is within the numerical relaxation region, in 
which the density and mass profiles are artificially shal- 
low. Thus, we tend to find fewer subhalos and satellites 
that satisfy the M300 and star-formation criteria. In the 



8 



50 
40 
30 
20 
10 
50 
40 
30 
20 
10 
50 
40 
30 
20 
10 




prq | | r- 

All nodes 

Dynamical friction 

^ M300 cut All 




I I I I I ll| — I M I lll'l 
Galaxies 




M300 Galaxies 



,V f 




1 10 100 1000 10" 

Number of Subhalos/Galaxies per Host 



FIG. 2: Distribution of the number of subhalos or satel- 
lites per halo in the host halo sample for CDM. The up- 
per panel shows the numbers of subhalos in the "all nodes" 
and "dynamical friction" samples both with and without the 
M300 > 5 x 10 6 Mq cut. The middle panel shows the distribu- 
tion in the number of subhalos satisfying the star-formation 
criterion for z rc = 7 and z IC = 11, and the bottom panel shows 
the distribution in the number of satellites satisfying both the 
M300 and star formation criteria. 



rightmost panels, the inner cut-off is set to r — 0.08i? V i r . 
Here we see that, because the mass profile is a bit steeper 
here than at r = 0.04i? v i r , we tend to overestimate M300, 
and hence we find that far more satellites and subhalos 
satisfy the criteria. Given that the mass profile ought to 
be becoming shallower inwards of r = 0.04i? V j r , we are 
often still overestimating M300, and hence the constraints 
on «k — t space based on the work in this section are quite 
conservative. 

In Fig. [5l we show our exclusion regions based on 
the number of star-containing satellites with M300 > 5 x 
1O 6 M0. In order to exclude a point in — r parameter 
space, we require that less than 5% of the host halos in 
the sample have at least 200 satellites that satisfy both 
the star-formation and M300 criteria. A point is allowed 
if at least 5% of the hosts have at least 200 satellites. 
One can see that all the CDM samples (Fig. [5]) have 
more than sufficient subhalos. 

We show constraints for the z rc = 7 (left) and z xe = 11 
samples (right). The region below the solid lines and 
to the right of the dashed line was previously excluded 
by the z — galaxy-cluster mass function and mass- 
concentration relation [30( ■ The Milky Way satellite lim- 
its are generally less constraining in that region. The 
light region to the left of the dashed line corresponds 



to constraints from the "all nodes" samples, and the 
dark red region corresponds to the additional region ex- 
cluded by the "dynamical friction" sample. Overlaid on 
both plots are the exclusion regions based on the sub- 
halo samples (i.e., without the star- formation criterion). 
The lower black line corresponds to the "all nodes" limit, 
while the upper line corresponds to the "dynamical fric- 
tion" sample. We find that star-formation criterion does 
affect the exclusion regions, although only at the level of 
a factor of ~ 2 in r for fixed unless < 20 km s" 1 . 

We also apply constraints from the highest-M3oo satel- 
lites, which are almost entirely the classical dwarf galax- 
ies for which the sample is currently complete. Of 
the 11 classical dwarfs, seven have M300 inferred to be 
M300 > 1O 7 M , one has a smaller M 30 o, while mass 
modeling is difficult for the remaining three and has so 
far precluded robust M300 estimates (Large and Small 
Magellanic Clouds and Sagittarius) [39|, |41|. Since the 
LMC and SMC have quite large v max , and there is a 
mild correlation of w max and M300, these two also likely 
have large M3oo- Thus, it is possible to constrain and 
t by determining the number of simulated satellites with 
M300 > 10 7 Mq, and making sure that the classical dwarfs 
are accounted for. When we perform this exercise, we 
find similar constraints as for the number of satellites 
with M 300 > 5 x 1O 6 M only if u k > 50 km s -1 . This is 
again because M300 is mildy correlated with w max , and in 
order for satellites with high v max to be affected strongly 
by decays, i>k needs to be significantly greater than v v - lr . 
However, for > 200km s _1 , constraints from the mass- 
concentration relation are typically stronger unless the 
z Ie = 11 "dynamical friction" sample characterizes the 
satellite population well. 



B. N(> Vmax) 

A standard way to characterize subhalos and satel- 
lites is by their velocity function, as illustrated in Fig. 
Q] [53|, Hij . This is often used in lieu of the subhalo or 
satellite mass because w max is relatively insensitive to the 
definition of the outer edge of the subhalo or satellite. 
Here, we explore the velocity function and the possibil- 
ity of constraints using the w max function of the observed 
satellite population. We emphasize that any constraints 
we find in this section are highly conservative because 
we determine w max for the subhalos and satellites in the 
absence of tidal stripping. Moreover, we compare the ve- 
locity functions in the decay parameter space with the ve- 
locity function of known dwarfs, corrected only for SDSS 
sky coverage (but NOT completeness). 

We consider a decay model to be ruled out if it fails to 
produce a sufficient number of satellites to reproduce the 
observed velocity function, but consider a model to be 
allowed if it overshoots the velocity function. In general, 
it is much easier to reduce i> max (e.g., by tidal stripping) 
than it is to increase it. 

We use Fig. [5] to illustrate a few salient and generic 
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FIG. 3: Distribution of the number of subhalos or satellites per halo in the host halo sample. Each plot shows distributions for 
fixed t>k and r (top row: Vk = 30 km s _1 ; bottom row: = 200 km s _1 ; left column: r = 20 Gyr; right column: r = 60 Gyr). 
Panels have the same meaning as in Fig. [2] 



features of the velocity functions. The velocity function 
of the "all nodes" subhalo population (denoted in Fig. [6] 
by the thick black line with 25% and 75% percentile bars) 
always lies above the velocity function for the observed 
dwarf galaxy population corrected for SDSS sky coverage 
(e.g., the dotted line in Fig. [1] which is that of Mada u 
et al. |1, compiled from data in Refe. [H ll0&Jl08| ). 
With a cut on the subhalo population of M300 > 5 x 
10 6 M Q (denoted by the thin black line with percentile 



bars), the velocity function constrains t ^ 20 Gyr for 
Vk ^ 200 km s _1 , which is only competitive with the 
decay constraint based on the number of subhalos above 
the M300 threshold for > 100 km s _1 . At such v^, 
it becomes nearly impossible to find subhalos the size of 
a Large or Small Magellanic Cloud in any Milky- Way 
mass halo because the decays greatly disturb the large 
subhalos that would have had large M300 in the absence 
of decays. For this range of v^, though, the observed 
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FIG. 4: Dependence of the subhalo/satellite distribution on the inner cut-off for the simulation interpolation table, for Vk = 
30 km s _1 and r = 20 Gyr. Left: Inner cut-off well within the numerical relaxation range. Center: Inner cut-off just outside 
the numerical relaxation range, as determined from simulations in which decay is turned off either at the beginning or later in 
the simulation. Right: Inner cut-off well outside the numerical relaxation region. Plot structure identical to that in Fig. [3] 
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FIG. 5: Exclusion limits in the — r parameter space. In both plots, the blue region marked "allowed" indicates the region of 
parameter space that has not yet been excluded. The red region marked "ruled out" and to the right of the dashed line is the 
part of parameter space that has been ruled out by observations of the galaxy-cluster mass function and the mass-concentration 
relation in galaxies, groups, and clusters [28-30] . The solid lines and the regions to the left of the dashed lines show limits 
from this work. The lower solid line shows the limit on the parameter space from the "all nodes" subhalo sample, and anything 
below the upper solid line is also excluded based on the "dynamical friction" subhalo sample. The light red region corresponds 
to limits using the sample of "all nodes" subhalos satisfying the star formation criterion with redshift z IC , while the dark red 
region shows the additional excluded region using the "dynamical friction" satellite sample satisfying the same star formation 
criterion. Left panel: z TC — 7 Right panel: z [c — 11. 



mass-concentration relation of galaxies rules out a greater 
swath of t. However, for t = 40 Gyr as illustrated in 
Fig. there is no problem forming sufficiently high- 
Vmax satellites. 

Since, in general, the velocity functions for both the 
z IC = 7 and z rc = 11 "all nodes" samples merge with the 
subhalo sample for large v ma x, the constraints on — r 
space are typically identical in the range in which con- 
straints from the velocity function are competitive with 
those found in Sec. IIII Al 



The velocity function of the "dynamical friction" sub- 
halo population produces constraints competitive with 
the constraints in Sec. IIII Al only for > 200 km s , 
° r ^ 100 km s _1 with the cut on M300, which we 
show in right-hand side of Fig. [6l For the subhalo and 
z IC = 7 satellite samples, this constraint again arises 
from the fact that it is difficult to produce the Magel- 
lanic Clouds. The constraint is tighter for the "dynamical 
friction" sample than for the "all nodes" sample because 
there are far fewer high-w max subhalos even in CDM due 
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to the fact that dynamical friction is more efficient for 
high-mass (and hence, high-w max ) subhalos. 

The velocity function provides the strongest constraint 
on r for > 30 km s _1 for the z ro = 11 "dynamical 
friction" satellite sample, such that r < 60 km s _1 is 
excluded for this «k range. The reason for this is ap- 
parent in Fig. |51 There is simply a dearth of satel- 
lites with u ma x > 10 km s" 1 . This limit is relatively 
insensitive to the M300 cut, since a large fraction of the 
z IC = 11 "dynamical friction" satellites have large M^oo- 
At fixed r, the velocity functions look quite similar for 
10 km s _1 < i>max < 20 km s~ 1 across a broad stretch 
of «k due to the fact that most of the subhalos are in 
the adiabatic regime of decay in which i>k 3> « V ir and 
t S> £dyn, where idyn is the typical dynamical time of par- 
ticles in the subhalos. We show the exclusion limits for 
z TC = 11, including both the constraints from N(> u ma x) 
and the results of Sec. MI Al in Fig. [SJ Note that these 
constraints are even stronger than for the galaxy-cluster 
mass function and the mass-concentration relation for 
Vk ^ 200 km s _1 , but only if the z rc = 11 "dynami- 
cal friction" sample is a good representation of the real 
satellite population. 

We note that stronger constraints are possible with 
N(> v max ) if one were to correct the observed velocity 
function for the SDSS selection function. 



C. M300 mass function 

Another possible way to constrain the decay parameter 
space is to use the full M 30 o mass function instead of the 
cuts we employed in Sec. IIII Al We show the CDM M300 
mass function for our subhalo and satellite samples in 
Fig. [5] We show the mass function for both "all nodes" 
and "dynamical friction" subhalo populations, as well as 
the corresponding satellite populations for z rc = 7 and 
z rc = 11. The cut-off in M 300 near M 300 = 10 6 M Q for the 
subhalo populations is an artifact of the mass resolution 
of the merger trees. There are several general features 
of this plot. First, most of the subhalos and satellites 
in our samples have M300 in the range corresponding to 
Milky Way satellites. Second, the location of the peak 
of the M300 mass function appears to depend somewhat 
on the star- formation prescription. The z rc = 11 satel- 
lite populations are skewed towards higher M300 than 
the z Te = 7 mass functions regardless of whether we con- 
sider the "all nodes" or "dynamical friction" samples. In 
addition, there are far fewer low-M 30 o satellites in the 
z IC = 11 samples than the z rc = 7 samples. Third, the 
high-M3oo tail depends quite strongly on whether or not 
dynamical friction is accounted for. There is a sharp 
cut-off in all the "dynamical friction" subhalo and satel- 
lite populations near M 30 o ~3x 1O 7 M . This cut-off is 
only slightly above the observed maximum M300 of the 
Milky Way satellites for which estimates of M300 exist. 
Moreover, the relatively narrow width of the M300 mass 
functions for the satellite populations suggests that the 



narrow range of M300 in observed satellites is a natural 
consequence of CDM cosmologies. The fact that the nar- 
row range of M300 in observed satellites can be simply 
explained in CDM has been previously noted by Stringer 
et al. [l09j . 

We compare the CDM M300 mass function to a few 
example M300 mass functions for decaying-dark-mattcr 
cosmologies, as shown in Fig. [§J We show mass functions 
for v k = 30, 100 km s^ 1 and r = 20,40 Gyr. The mass 
functions can look quite different, depending on the decay 
parameters. In general, smaller lifetimes lead to a much 
broader smearing of the mass function, with a low-M 30 o 
tail predicted for satellite as well as subhalo populations. 
While the peak of the mass function necessarily shifts 
to lower M300 as r gets smaller, it is less sensitive to 
variations in v^, although the shape of the mass function 
clearly is quite sensitive. As Wk increases, the high-M3oo is 
more sharply cut off since even the largest subhalos begin 
to become highly disturbed as a result of the decays. 

While these plots show that M 30 o mass functions de- 
pend on decay parameters, they also show that they de- 
pend on the star-formation prescription and on dynam- 
ical processes once the subhalo or satellite is accreted 
onto a halo. This is apparent in Fig. [§1 and what is es- 
pecially striking is the high-Msoo tail of the "all nodes" 
samples which extend to far higher M300 than the "dy- 
namical friction" or observed satellite samples. In addi- 
tion, one effect we have not modeled is tidal stripping. 
While 300 pc is typically deep within a satellite or sub- 
halo, the mass within that radius may be reduced as a 
consequence of tidal stripping, as high-apocenter orbits 
are progressively removed from the center. However, if 
M300 were to be significantly affected by tidal stripping, 
the associated stellar population should also look signif- 
icantly disturbed, unless the dark-matter orbits are, on 
average, highly radial. Although we hypothesize that 
the "dynamical friction" samples are more likely to rep- 
resent reality than the "all nodes" samples, the shape of 
the M300 mass function will still depend on the physics 
of star formation and tidal stripping. The latter can in 
principle be modeled using cosmological iV-body simula- 
tions, but would require extremely high resolution. 

There are a few things that are likely to be robust 
to these effects. The first is the upper M300 tail of the 
mass function, since while it is possible to remove mass 
from inside 300 pc in a halo, it is hard to add mass if 
there are also only negligible amounts of baryons in the 
halo (which might have compressed the dark-matter mass 
profile) . This is essentially what we discussed at the end 
of Sec. IIII Al in which we used the classical satellites 
(which tend to have M 30 o ~ 10 7 M Q ) In Fig. U3 we see 
that for «k = 100 km s _1 and r = 40 Gyr, we expect no 
satellites with M 300 >2x 10 7 M Q ; if ever such a dense 
satellite were discovered, it would rule out that point in 
decay parameter space. 

Thus, some aspects of using the M300 mass function to 
constrain the nature of dark matter are more tractable 
than others. A sky- and selection function-corrected 
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v k = 100 km/s, t= 40 Gyr, all nodes v k = 100 km/s, t= 40 Gyr, dynamical friction 




FIG. 6: Maximum circular velocity functions for i>k = 100 km s^ 1 and r = 40 Gyr, for "all nodes" subsamples (left) and 
"dynamical friction" subsamples (right). The line types indicate different subsamples as indicated in the legend: all subhalos, 
subhalos satisfying the star formation criterion with z la = 7, and subhalos satisfying the star formation criterion with z la = 11. 
The thin lines of each type indicate that a cut of M300 > 5 x 10 6 Mq has been included, and thick lines show the « ma x distribution 
without a cut on M3oo- 



z re = 11 satellite + large-halo limits 




FIG. 7: Exclusion limits for the z la = 11 satellite populations 
including the limits from the velocity function. Lines and 
shading of the plot have the same meaning as in Fig. [5] 



M300 distribution function has not yet been published, 
but there is at least one group working on this (Wolf et 
al., in prep.). 
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FIG. 8: M300 distribution for CDM. Line types have the same 
meaning as in Fig. [6] Thick lines represent "all nodes" sub- 
samples, and thin lines represent "dynamical friction" sub- 
samples. 
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FIG. 9: M300 distributions for those consistent with limits on the t parameter space using the "all nodes" subhalo sample 
(top) and those consistent with the "dynamical friction" subhalo sample (bottom). The line types have the same meaning as 
in Fig. [6] The thick lines represent "all nodes" and thin lines represent "dynamical friction" subsamples. 



IV. DISCUSSION 

In this work, we have shown that the low-i>k end of the 
Uk— t decay parameter space is currently best constrained 
by the number of satellites with M300 > 5 x 10 6 M Q . 
We have found that the precise constraint depends on 
the star-formation prescription, but not dramatically so 
(within a factor of 3 for r for fixed v^). We found that 
the velocity function of satellites and the M300 values 



of the classical dwarfs provide similar constraints for 
vy_ > 100 km s _1 , which is because decays only affect 
the largest subhalos for such v^. However, in the case 
of the 2 rc = 11 "dynamical friction" sample, the velocity 
function more strongly constrains the decay space than 
the number of satellites above the M 30 o threshold, the 
mass-concentration relation, or the galaxy-cluster mass 
function. The limits we set in Figs. [5] and [7] are quite 
conservative because we use merger trees with high erg 
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and we use the subhalo properties at accretion to deter- 
mine Wmax and M300 at z = 0. 

We showed that the distribution of M 30 o might be an 
avenue for future better constraints of decay parameter 
space, although the distribution does appear to depend 
on the details of star formation in subhalos. Conversely, 
one can think of the M300 distribution as a way to probe 
star-formation physics as well as dark-matter physics. 
In addition, tidal stripping can potentially lower M 30 o, 
although typically 300 pc is smaller than the radius at 
which the circular velocity curve of subhalos peaks, mak- 
ing M300 likely a more constant property of a subhalo 
over its lifetime than w ma x- Moreover, if M300 were to be 
significantly affected by tidal stripping, there would likely 
be evidence of tidal stripping in the stars, too. There are 
some generic features in the M300 distributions that likely 
to be robust, such as the high-M 30 o tail of the distribu- 
tion. 

In theory, M300 provides a cleaner and potentially more 
powerful probe of dark-matter properties than w max be- 
cause !) max typically occurs at larger radii than 300 pc, 
at least for CDM cosmologies. The fact that M300 
probes the innermost mass of the satellite galaxies is in- 
teresting for dark-matter theories that are alternatives 
to CDM. All non-CDM theories invoke energy injection 
or transfer i nto t he dark-matter population by decays 
(e.g., flol. [24l [TToj ) or by introducing a non-trivial colli- 
sion term into the dark- matter Boltzmann equation (e.g., 
[ill llllj ). Typically, the effects of such energy injec- 
tion or collision has been parametrized by Q, the coarse- 
grained dark-matter distribution function. Q is enor- 
mous for CDM, but becomes is modest once decays or 
collisions are turned on. 

The problem with trying to infer Q from data is that Q 
depends on the velocity structure of dark matter, which 
can never be directly measured (although it may be in- 
directly measured if, for example, the annihilation cross 
section is velocity-dependent). Specific non-CDM the- 
ories predict relations between Q and the dark-matter 
density profile, but one must analyze the data in the 
context of that specific model [l(J Hj] . While this is use- 
ful to constrain specific theories, one cannot generically 
determine if the observations deviate CDM on the basis 
of Q. 



The advantage of using M300 to characterize the dark- 
matter halos, and using that to consider deviations from 
CDM, is that any non-CDM dark-matter theory implies 
that M300 is lower than the CDM value. Due to the nega- 
tive heat capacity of self-gravitating halos, kinetic energy 
injection or collision terms tend to reduce the central den- 
sity of the halos (at least until core-collapse, in the case 
of collision terms). Thus, any physics that would reduce 
Q would also reduce M 30 o- Since 300 pc is typically deep 
within the halo, it could be a good probe of dark-matter 
properties. 

The inferred property of the dark content of Milky Way 
satellites that has the smallest errors is the m ass w ithin 
the half-light radius of the stars, M 1/2 0, [H, OH • The 
half-light radius is often smaller than 300 pc (e.g., ~ 
100 pc for Coma Berenices [HI, [Hj]), although it is often 
larger, especially for the classical dwarfs. However, the 
small value of the half-light radius for some of the dwarfs 
is actually beneficial in constraining the nature of dark 
matter, since one is measuring the mass within a tiny 
radius centered on the potential minimum of the halo. 

If the M300 mass function is to be a useful probe of ei- 
ther dark-matter physics or galaxy evolution, it is neces- 
sary to see how much M 30 o is affected by the dynamics of 
subhalos inside the host halos. However, if this is possible 
to determine, then in principle, the M300 mass function 
or the M1/2 mass function of satellites should provide an 
interesting window into dark-matter properties. Upcom- 
ing wide-field surveys, such as the Dark Energy Survey, 
SkyMapper, and the Large Synoptic Survey Telescope, 
should reveal many more Milky Way satelli tes and m uch 
more about dark-matter properties |U EH, Ill3l - lll5j | . 
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